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Abstract We show the exact solution of Bogoliubov equations at zero-energy in the 
critical supercurrent state for arbitrary shape of potential barrier. With use of this 
solution, we prove the absence of perfect transmission of excitations in the low-energy 
limit by giving the explicit expression of transmission coefficient. The origin of disap- 
pearance of perfect transmission is the emergence of zero-energy density fluctuation 
near the potential barrier. 
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1 Introduction 

In 2001-2003, Kovrizhin and his collaborators^ have shown an interesting theoretical 
prediction on the tunneling properties of Bogoliubov excitations [2]; the Bogoliubov ex- 
citations show perfect transmission across a potential barrier in the low-energy limit. 
It is called anomalous tunneling. It is quite different from an ordinary particle obeying 
Schrodinger equation, which shows perfect reflection in the low-energy limit. Later, 
Danshita et al. [3] have extended the problem in the presence of the supercurrent. Solv- 
ing the delta-functional barrier problem, they have shown that (a) perfect transmission 
occurs even when the condensate supercurrent exists, except for the critical supercur- 
rent state; (b) under the critical supercurrent, the partial transmission occurs. Thus, a 
consistent explanation for both perfect transmission in the non-critical states and the 
absence of perfect transmission in the critical state had been highly desired. As for the 
perfect transmission in the non-critical states, its physical mechanism has been investi- 
gated in many works 4 , 5 and importance of the similarity between the condensate and 
low-energy excitations has been pointed out. However, few results have been known on 
the critical supercurrent state. It is crucial to understand the reason why the critical 
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supercurrent state shows an exceptional behavior, for the purpose of constructing a 
unified picture for non-critical and critical states. In our study, we prove the absence 
of perfect transmission in the critical supercurrent state for arbitrary shape of barrier, 
and clarify its physical mechanism. 



2 Formulation 

We begin with time-dependent Gross-Pitaevskii(GP) equation 

iJ^(M) = (-5^2 + tW) V'OM) + \il>(x,t)\ 2 1>(x,t). (1) 
Here we use a dimensionless description. Assuming the solution in the form of 

i>(x,t) = e~ iflt ^(x) + [«(i)e- Ut -/(i)c id ]}, (2) 
and ignoring the higher-order terms of u and v, we obtain stationary GP equation 

W{x) = 0, L = -i^ + u ( x ) - P + l 1 ^)! 2 ( 3 ) 
for the condensate wavefunction, and Bogoliubov equations 



u(x)\ _ ( u(x) 



-(&(x)*) 2 L+ \&(x)\ 2 J \v(x)J \-v(x) 



(4) 



for wavefunctions of Bogoliubov excitations. 

Henceforth, we would like to consider the problem depicted in Fig. [1] Therefore, 
we assume that U(x) is a short-ranged potential barrier, and that the condensate 
wavefunction has the following asymptotic form: 

ip(x — » ±00) = exp ^i^qx ± ^ + const. . (5) 

This asymptotic behavior determines chemical potential as /i = 1 + q 2 /2. Setting 
'I'(x) = A(x) exp[i0(x)], one obtains 

1 d 2 , TT , q 2 ( 1 ,\ , , ,2 



HA = 0, H = --— 2+ U + ^[ J ,-l)-l + A^ (6) 

Here the second equation is already integrated once, and a constant of integration 
becomes q from Eq. (J^}. Figure [2] shows the Josephson relation. From this figure, the 
condition for the critical supercurrent state is given by 

!=°- < s > 

By means of the condensate phase 0, we introduce the following quantities: 

ib = we + ve , G = tie — ve . (9) 



3 



*(£-►— oo) = 



n exp[i(gx + |) 




Fig. 1 Supercurrent state of the condensate wave- 
function and tunneling problem of Bogoliubov exci- 
tations across a short-ranged potential barrier. The 
shaded area shows the condensate density profile 
\ty(x)\ 2 , and q represents the magnitude of the su- 
percurrent (not a wavenumber of a Bogoliubov ex- 
citation). Such a reflectionless supercurrent solution 
must have the phase difference ip, and supercurrent q 
depends on tp as shown in Fig. f2] In this supercurrent 
state, we further consider the tunneling problem of 
Bogoliubov excitations, i.e., find the transmission 
and reflection amplitudes t and r for an incident 
wave p^l x 




Fig. 2 Josephson relation q(ip) for a 
delta functional barrier U(x) = Vb<5(x) 
with V = 1, 5, and 10. Though the 
exact expression of q(tp) is compli- 
cated, when Vb 2> 1, it approximately 
becomes q(tp) ~ sin tp + 0(^^) [3J. 

(ip c ,q c ) are the critical supercurrent 
states, where the supercurrent reaches 
to the maximum value. An example 
of Josephson relation for a rectangular 
barrier can be seen in Ref. [6]. 



Bogoliubov equations are then rewritten as 



* 8 -5=(f)-« 

f + "*)°-5=(3)-« a 



(10) 
(11) 



S and G can be interpreted as phase and density fluctuations, because one can show 
the following expressions from Eq. ([2]): 

2 



A 2 [l + - Re(Ge- iet )] + 0(S, G) 2 



= e -^ t+ie [l + 1 Im(Se- iet )] + 0(S, Gf 



(12) 
(13) 



In order to obtain the transmission amplitude, we construct the tunneling solution 
with the following asymptotic form: 



1 

fc?/[2(e - qk x )) 
1 



ifci x . 
C + 



te 



1 



fc|/[2(e - qk 2 )] 



re (x —> —oo) 



(14) 



(x — » +co). 



Here k\ and k<i are real positive and negative roots of the dispersion relation e 
qk+\sfk 2 (k 2 +4). 



3 Exact Solution of Bogoliubov Equations at e = in the Critical 
Supercurrent State for Arbitrary Shape of Potential Barrier 

Since Bogoliubov equations are two-component and second-order linear differential 
equations, the general solution can be expressed in terms of four linearly-independent 
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Fig. 3 Localized density fluctuation solution G = dA/d(fi for a delta functional barrier 
U(x) = Vq5(x) with Vq = 4.6. (We note that dA/dtp in the figure is multiplied by a constant.) 



solutions. We have obtained the exact solutions of Bogoliubov equations at e = valid 
only for the critical current state, in other words, only when the condition ((HJ holds. 
They are given by [8] 



-2iqAA : 

Au, 



(Sn 



(Sxv 



' ^ „• A f 

-2iqA v 



G u dx \ 
A* > 



Azdx 



) 



(15) 



V 



A, 



o A* 
da; 



o Al J 



Here we have introduced the following notations: 



OA 

dip ' 



A 3 (x) := 



Aip(x')&x' 
A(x' 



A3 



(16) 



(S\,G\), which can be expressed as (u,v) = (•J', is a well-known solution [7]. The 
solution (Sin, Gin), which can be expressed as (u,v) = (d'I'/dip, — d^* /dip), is specific 
to the critical current state, and represents the localized density fluctuation near the 
potential barrier. See Fig. [3] (Su,Gu) and (Siy , G\y) are exponentially divergent 
solutions. However, as we will see in the next section, all four solutions are necessary 
to prove the absence of perfect transmission. 



4 Absence of Perfect Transmission: Sketch of Proof 

Using the exact zero-energy solution shown in the previous section, we can prove the 
absence of perfect transmission. In the following, we show the sketch of proof. The 
detailed calculation is given in Ref. [8]. Henceforth, we assume U(x) = U(—x) for 
simplicity. 

i) Our goal is to obtain the tunneling solution (|14p up to first-order in e, that is, 

f 1 + - (0) + e (r~« + + x) + 0{<?) (x -+ -oo) 

S(x) — > I (17) 

[ i (0) +e (i(l) + _gi x ) +O(e 2 ) (X^+OO). 

Herei = i (0) +ei (1) + --- , and f = f (0) +ef (1) + ■ ■ ■ . 

ii) In order to achieve the above purpose, we construct the solution of Bogoliubov 
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equations in the form of power series with respect to e: (S, G) = XnLo e n (S^ n \ G^ r 
We then obtain the following inhomogeneous differential equations: 

-,(«) _ n d / G (n) \ _ „(„-!) 



H * n '-tt\TA-) = ^ l '> (18) 

(H + 2A 2 )G^-'^(^P i = S^- 1 \ (19) 

A CLX \ i\ J 

where the right hand sides should read as zero if n = 0. Since four homogeneous solu- 
tions are already given in Eq. (|15p . a particular solution can be found by the method of 
variation of parameters. Thus, (S^ n \ G*-™') can be determined from 
recursively. 

iii) Bogoliubov equations with finite energy e have four linearly-independent solu- 
tions. Generally, two of the four behave as plane waves far from the potential barrier, 
and the other two are unphysical solutions which diverge exponentially. In solving the 
tunneling problem, we are particularly interested in the former two solutions. 

iv) Therefore, we extend the non-divergent zero-energy solutions, i.e., (Si,G{) 
and (Sin, Gin), to first-order in e. The particular solution of Eqs. (|18p and fj 19p ob- 
tained by the method of variation of parameters diverges exponentially, so we must 
cancel the divergent term by adding the divergent homogeneous solutions (Sh,Gh) 
and (Siv,Giv). This manipulation is most important. After the calculation, we can 
obtain the asymptotic forms of the solutions extended up to first-order as follows: 

S\° tal (x) -^l + e( . q ~ T, x + jsgax) + 0(e 2 ), (20) 
Vig(l-r) ) 

const, x S|n al (z) — > sgnz + e f -^fT^sy kl + A j + Q(e 2 ). (21) 



Here 7 and A are constants, and r\ is defined by the following integral: 

1 OA 



V 



dxA— 

dip 



U0 



dx A3 dip 



(22) 



v) By making a linear combination of 5'* otal and 5'}'f I tal , we can construct the 
solution of the form (|17[) . The transmission amplitude is obtained explicitly as 

^=*( e -0) = -^_. (23) 

Unless r\ 7^ ±g, < l^ -*] 2 < 1 holds. Thus, the absence of perfect transmission is 
proved. We have confirmed in the delta-functional barrier model that r\ = q occurs 
only when there is no potential barrier. 



5 Discussion 

Let us consider the physical origin of the absence of perfect transmission. The tunneling 
solution in the low-energy limit can be written as [8] 
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This is drastically different from the non-critical states, in which the wavefunctions of 
excitations coincide with the condensate wavefunction[4][9]. In the non-critical states, 
there is only one non-divergent solution (u,v) = (if - , ]P*) <-> (S,G) = (A,0) at e = 0, 
so G cannot contribute to the wavefunction of excitations in the low-energy limit. 
In the critical supercurrent state, on the other hand, another non-divergent solution 
(u, v) = (d\P/dp, — d^P* /dp) <-> (S,G) = (~2iqAA^, Atp) arises, and contributes to the 
low-energy wavefunction of excitations. As shown in Fig. [3] this solution represents 
the density fluctuation localized near the potential barrier. Thus, the presence of local 
density fluctuation near the barrier in the low-energy limit is the origin of the absence 
of perfect transmission. 

We further note that the critical supercurrent states are at the "phase boundary" 
which separates steady flow states and nonstationary flow states. (See, e.g., Fig. 1 
in Ref. [ID]-) Therefore, our result suggests that the emergence of low-energy density 
fluctuation can characterize the destabilization of the superflow. In Ref. [IT], it has 
been shown by numerical simulation of time-dependent GP equation that if q > q c , 
soliton-phonon creation occurs and stationary superflow no longer exists. We expect 
that the local density fluctuation G = OA/ dip possesses the information on the "manner 
of collapse of stationary superflow" . That is, the profile of the density fluctuation could 
describe a part of collapse phenomena, such as the soliton-phonon creation. Thus, the 
investigation of the role of low-energy density fluctuation near the critical supercurrent 
state is left as a future work. 

6 Conclusion 

In conclusion, with use of the exact solution of Bogoliubov equations at zero-energy, we 
have exactly proved the absence of perfect transmission of Bogoliubov excitations in 
the critical supercurrent state for arbitrary shape of potential barrier. The origin of the 
absence of perfect transmission is the emergence of low-energy density fluctuation near 
the potential barrier. Because of this density fluctuation, wavefunctions of excitations 
does not coincide with the condensate wavefunction in the low-energy limit. 
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